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Abstract 

We examine theoretically coherent electron transport through the single- molecule magnet Mni2 , 
bridged between Au(lll) electrodes, using the non-equilibrium Green's function method and the 
density-functional theory. We analyze the effects of bonding type, molecular orientation, and 
geometry relaxation on the electronic properties and charge and spin transport across the single- 
molecule junction. We consider nine interface geometries leading to five bonding mechanisms and 
two molecular orientations: (i) Au-C bonding, (ii) Au-Au bonding, (iii) Au-S bonding, (iv) Au-H 
bonding, and (v) physisorption via van der Waals forces. The two molecular orientations of Mni2 
correspond to the magnetic easy axis of the molecule aligned perpendicular [hereafter denoted as 
orientation (1)] or parallel [orientation (2)] to the direction of electron transport. We find that 
the electron transport is carried by the lowest unoccupied molecular orbital (LUMO) level in all 
the cases that we have simulated. Relaxation of the junction geometries mainly shifts the relevant 
occupied molecular levels toward the Fermi energy as well as slightly reduces the broadening of 
the LUMO level. As a result, the current slightly decreases at low bias voltage. Our calculations 
also show that placing the molecule in the orientation (1) broadens the LUMO level much more 
than in the orientation (2), due to the internal structure of the Mni2. Consequently, junctions 
with the former orientation yield a higher current than those with the latter. Among all of the 
bonding types considered, the Au-C bonding gives rise to the highest current (about one order 
of magnitude higher than the Au-S bonding), for a given distance between the electrodes. The 
current through the junction with other bonding types decreases in the order of Au-Au, Au-S, 
and Au-H. Importantly, the spin-filtering effect in all the nine geometries stays robust and their 
ratios of the majority-spin to the minority-spin transmission coefficients are in the range of 10^ to 
10^. The general trend in transport among the different bonding types and molecular orientations 
obtained from this study may be applied to other single-molecular magnets. 

PACS numbers: 73.63.-b, 85.65.+h, 85.75.-d, 75.50.Xx 
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I. INTRODUCTION 



Advances in experimental techniques have led to a great number of experimental stud- 
ies on electron transport through molecular junctions formed by single molecules bridged 
between electrodes or molecular monolayers adsorbed onto surfaces, using three-terminal 
set-ups or scanning tunneling microscopy (STM) measurements. Recently, molecular junc- 
tions based on single-molecule magnets (SMMs) connected to electrodes or monolayers of 
SMMs at surfaces, have been fabricated and their electron transport characteristics^"- have 
been measured, as well as their mechanical, electronic, and magnetic properties^i"— . Elec- 
tron transport through an SMM drew a lot of attention because of the intriguing interplay 
between its transport properties and the internal magnetic degrees of freedom, which is ab- 
sent in transport through small organic molecules. An SMM consists of several transition 
metal ions interacting through organic or inorganic ligands via super-exchange interactions. 
The spin configuration of the ground state of an SMM is determined by a delicate balance 
among the super-exchange interactions of different strengths between the transition metal 
ions. Thus, the magnetic structure of an SMM must be taken into account in understanding 
its electron transport and other quantum properties. Recently, first-principles calculations 
of transport through an SMM were performed on a prototype SMM Mni2 terminated with 
a thiol (-S) group within Au electrodes.—"— The calculations carried out in Refs. 1^. sug- 
gest that the Mni2 molecule can function as a spin filter with low bias voltage even with 
non-magnetic electrodes. Even though the Mni2 is chemically bonded to the Au electrodes, 
the broadening of the relevant molecular orbitals is so small compared to its charging en- 
ergy that the Kondo temperature is expected to be extremely low. In addition, these^^i^ 
and other calculationsi^ii^ demonstrate the importance of the internal magnetic degrees of 
freedom of the Mni2 in electron transport, in contrast to typical quantum dots, regardless 
of the specific details of the coupling of the molecule to the electrodes. 

Systematic studies on molecular junctions based on small organic molecules reveal that 
the current through the molecular junctions is highly sensitive to properties of interfaces 
between the molecules and electrodes because the interfaces determine the degree of the 
overlap between the molecular levels and conduction channels of the electrodes.—"— A thiol 
group is most commonly used to build a strong chemical link between the single molecules 
and the Au electrodes in molecular junctions. The conductance of such single-molecule 
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junctions is typically several orders of magnitude smaller than Go = ^2.6^ /h (the conductance 
quantum). Binding of the S atoms to the hollow sites of the Au surface gives rise to differ- 
ent conductance (by a factor of about 2 to 3) from binding of the same S atoms to on-top 
sites of the Au surface.— Molecular junctions terminated with an amine (-N) group within 
Au electrodes reveal even lower conductance than those with a thiol group. In order to 
increase conductance of molecular junctions, linker molecules including Au atoms^^ or iso- 
cyanide derivatives^^ were used between the single molecules and electrodes. In addition, 
recent experiments show that conductance through molecular junctions based on small sin- 
gle molecules can be enhanced to an order of Gq by their direct bonding to the electrodes 
without a thiol group or any other linker molecules.— This enhancement is attributed 
to a strong coupling between the molecules and the electrodes, which places the transport 
in the transparent regime rather than in the tunneling regime.— However, corresponding 
systematic studies have not yet been carried out for electron transport through an SMM. 
In the case of quantum dots, the properties of interfaces are negligible in transport. How- 
ever, for molecular junctions based on SMMs, the molecules are chemically bonded to the 
electrodes, and thus their transport properties can change with different interfaces even if 
the molecules are weakly coupled to the electrodes in the sense that the charging energy is 
much greater than the broadening of the relevant molecular levels. In this paper, we inves- 
tigate how bonding types, interface geometries, and geometry relaxation influence transport 
characteristics of an SMM Mni2 when it is bridged between Au electrodes. 

An SMM Mnia^O consists of four inner Mn^"*" ions {S = 3/2) surrounded by eight outer 
Mn^"*" ions {S = 2) which are all antiferromagnetically coupled through O anions, as well 
as C and H atoms, as shown in Fig. [H The dominant exchange interactions are between 
the inner Mn^+ ions and the outer Mn^"*" ions.™ In the ground state the magnetic moments 
of the four Mn^"*" ions are antiparallel to those of the eight Mn^"*" ions, such that the total 
spin becomes S' = 2x8-3/2x4 = lO.-^"-^ An SMM Mnu initially synthesized by Lis^ 
consists of ligands which do not form direct chemical bonding to an Au surface. Thus, most 
transport experiments through an Mni2 within Au electrodes were performed on systems 
where the Mni2 molecules are chemically bonded to surfaces or electrodes via linker molecules 
such as a thiol group,-"^!^ or physically adsorbed onto surfaces or electrodes without any 
linker molecules.- To form direct Au-S bonding, some of the original ligands of the Mni2 
molecules^ were substituted by S-containing ligands, or the Mni2 molecules were deposited 
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FIG. 1: (Color online) Top view of an isolated Mni2 molecule: Mn (purple), O (red), C (gray), H 
(white). The ground-state spin is = 10 or the magnetic moment is 20 /i^. 

onto an Au surface that was initially functionalized with S-containing alkane chains.- So 
far, the measured electric current through an Mni2 molecule or its derivative in various 
experiments is in the range of 1 pA to 100 pA at bias voltage of a few tens of mV.-i^ 
Compared to binding a small organic molecule to Au electrodes through a thiol or amine 
group, binding an Mni2 to Au electrodes via linker molecules bears the following differences: 
(a) an Mni2 has much more binding sites to the linker molecules, and (b) the orientation of 
an Mni2 relative to the electrodes is, to a great extent, determined by binding sites of the 
linker molecules to the Mni2, rather than their binding sites to the electrodes. 

Full control of the properties of interfaces is extremely difficult to achieve in molecu- 
lar junctions and tunneling measurements using STM or scanning tunneling spectroscopy 
(STS).— i^^i""— i^liSS Thus, in transport experiments, for a given bonding type, several binding 
sites are plausible, and hundreds of fabricated samples of the bonding type give rise to a 
histogram of conductance for given gate voltage or a series of current-voltage curves. In the 
present study, we take into account five different bonding types between an Mni2 and Au 
electrodes, and for each bonding type, some representative interface geometries are exam- 
ined. Even though we do not simulate all possible configurations of the interface that may 
be realized in fabricated samples for a given bonding type, the general trend in our calcu- 
lated transport properties will elucidate the effects of linker group, molecular orientation, 
and geometry relaxation, as well as the effectiveness of different bonding types and interface 
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geometries in transport through an Mni2 and other SMMs. 

In this work, we consider an Mni2 molecule bridged between Au(lll) electrodes (two- 
terminal set-up) via nine different ways (Fig. and investigate their current- voltage char- 
acteristics using the density-functional theory (DFT) and the nonequilibrium Green's func- 
tion method (NEGF). The nine different ways can be categorized according to their bonding 
types: (i) Au-S bonding, (ii) Au-C bonding, (iii) Au-Au bonding, (iv) Au-H bonding, and 
(v) physisorption through van der Waals forces. DFT does not fully capture the nature of 
van der Waals forces where the orbitals considered do not overlap. We note, however, that if 
molecular orbitals do not overlap with the orbitals of electrodes at a molecular junction, the 
transport properties of the junction belong to the tunneling regime. A recent study reveals 
that first-principles calculations based on DFT and NEGF provide qualitative features of 
transport in such a regime.— For example, when a Pt break junction becomes broken with 
vacuum between the electrodes, a first-principles calculation showed exponentially decaying 
conductance as a function of distance between the electrodes.— This result suggests that our 
methodology can be used to understand the transport in the case of physisorption, although 
the physical mechanism in our case differs from that in the broken junction. In the present 
calculations, we consider two molecular orientations relative to the electrodes: when the 
magnetic easy axis of the Mni2 is normal [orientation (1)] or parallel [orientation (2)] to the 
transport direction, the z axis, as shown in Fig. |2l The structure of the Mni2 dictates that 
in the orientation (1) one linker molecule can be attached to each electrode without signif- 
icantly deforming the Mni2, while in the orientation (2), two linker molecules are possibly 
accommodated for each electrode without much distortion. As illustrated in Fig. |2l the nine 
different interface geo metries are named based on the following rule. Adapted from the nota- 



tions used in Refs. 



14 



15|, we refer to the orientation (1) as Geo 1, and to the orientation (2) 



as Geo 2. For each interface geometry, the orientation of an Mni2 is first specified, then the 
type of adsorption (either chemical or physical) is given. For chemical bonding, the type of 
bonding and the name of the linker group are noted, while for physisorption, whether linker 
molecules are placed or not is stated. For example, Geo l:Au-(SC3H6)2-hollow [Fig. EJ^a)] 
denotes an interface geometry with the Mni2 in the orientation (1) and the linker molecule 
(SC3H6)2 bonded to the hollow sites of the electrodes via Au-S bonding. We will first out- 
line our model and computational method. Then we will discuss the interface geometries 
in detail, and the effects of geometry relaxation, molecular orientation, bonding type, and 
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linker group on transport properties of the Mni2. A brief conclusion will follow. 



II. COMPUTATIONAL METHOD AND MODEL 

We use the quantum transport code SMEAGOL^'^ interfaced with the DFT code SIESTA^S. 
The Au electrodes are treated semi- infinitely using SIESTA. The scattering region consists 
of the Mni2, linker molecules, and six flat Au atomic layers on each side of the Mni2 (6x6 
surface atoms per layer), as shown in Fig. |2J Even if the number of the Au atomic layers 
increases from six to nine on each side of the Mni2, the transmission probability does not 
change at all. This was tested and confirmed for the interface geometry with the Au-C bond- 
ing [Geo l:Au-(C3H6)2, Fig. |2]^c)], where the broadening of the relevant molecular level turns 
out to be the largest among the interface geometries considered. The electronic structures of 
the electrodes and scattering region are calculated within the spin-polarized Perdew-Burke- 
Ernzerhopf (PBE) generalized-gradient approximation (GGA)^, using SIESTA. We generate 
Troullier-Martins pseudopotentials^ for Au, Mn, S, O, C, and H with scalar relativistic terms 
and core corrections except for H. We also build corresponding basis sets of the elements 



using the scheme presented in Ref. 4l|. For the Mn basis set, 3p orbitals as a semi core 
are required in order for the molecular orbital levels of an isolated Mni2 to be comparable 
to those obtained from all-electron calculations^^. With these pseudopotentials and basis 
sets, self-consistent DFT calculations including spin-orbit coupling^ are performed for an 
isolated Mni2. Using a modified version^ of SIESTA, we obtain the total magnetic moment 
of 20 fiB and the magnetic anisotropy barrier of 66.4 K,^^. This is in good agreement with 
experiment^ and with the barrier— computed using the DFT code VASP^. To reduce the 
computational cost, a small Au basis set of a single s orbital is used for transport calcula- 
tions, while a large Au basis set of both d and s orbitals is used for geometry relaxation. 
The transmission at the Fermi level in Au is dominated by the s states because the d orbital 
levels are well below the Fermi level. It is checked that in the scattering region, the density 
of states projected onto Mn d orbitals with the large Au basis set is fairly similar to that 
with the small Au basis set, which justifies the utilization of the small Au basis set in the 
transport study. For a given interface geometry, the distance d between the electrodes was 
determined such that all of the terminating atoms from the linker molecules (such as S, C, 
or Au atoms) are bonded to the lowest-energy sites of the the Au slab with reasonable bond 
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FIG. 2: (Color online) The nine geometries of the scattering region with different bonding types 
and linker groups. Not all of the Au atomic layers are shown, (a), (b), (e), (f): Au-S bonding 
with different linker molecules and in the orientations (1) and (2). (c): Au-C bonding, (d): Au-Au 
bonding, (g) and (i): physisorption with and without linker molecules, respectively, (h): Au-H 
bonding. On the top and in the middle, the positions of the twelve Mn ions in relation to the 
binding sites for the linker molecules are labeled for the orientations (1) and (2), respectively. 
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lengths. For example, for the Au-S bonding, the distance d was obtained from an optimum 
bond length between a S-terminated alkane chain and the hollow site of a small Au clus- 
ter. In constructing the initial geometry of the scattering region, each component of the 
region, such as the electrodes, the Mni2, and the linker molecules, are separately optimized 
in advance. The largest force is found at the interface between the linker molecules and 
the electrodes.^^ Thus, the geometry relaxation with a fixed d allows the linker molecules 
and the Mni2 to relax further and lowers the force at the interface. Prior to the transport 
calculations, all of the geometries of the scattering region or interface geometries considered 
are relaxed with a fixed distance d between the electrodes, using SIESTA, until the magni- 
tude of the maximum force exerted on the atoms becomes less than 0.1 eV/A, unless stated 
otherwise. 

To avoid any quantum confinement effects,"2S for the scattering region, periodic bound- 
ary conditions are applied in the transverse directions and 3x3x1 fc-points are sam- 
pled. The retarded Green function for the scattering region, G^^, has the form of 
[e+^EM — -f^EM — Sl — where Sem and Hem are the overlap and Hamiltonian matrices 

for the scattering region, respectively.— Here Sl and S^; are self-energies arising from the 
interactions of the Mni2 with the left and right electrodes, respectively. These self-energies 
depend on energy E and k vector, and they are obtained from spin-polarized calculations 
of retarded surface Green functions for the electrodes which are computed using the scheme 



developed in Ref. 46|]. The initial magnetic moment of the scattering region is set to 20 
fiB- Interactions with phonons or additional electron correlations such as on-site Coulomb 
repulsion (Hubbard-like U term) are not taken into account. The density matrix of the 
scattering region is self-consistently solved within the NEGF formalism until it converges. 
After the convergence of the density matrix, the total magnetic moment slightly increases, 
such as 20.3 fiB for Geo l:Au-(SC3H6)2-hollow, due to a small amount of spin polarization 
in the Au atomic layers caused by the Mni2. Then the transmission T{E,Vb) at low bias 
voltage Vb is calculated as follows: 

T{E, Vb) = Tt[Tl Ti? G^u]iE, H), (1) 

where Tl and Tb denote the broadening of molecular levels induced by coupling to the left 
and right electrodes. Accurate calculations of T{E, Vb) require high resolution in energy E 
(such as 7/ieV for very sharp transmission peaks) due to the weak coupling between the Mni2 
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and the electrodes. Since the electrodes are treated semi-infinitely and periodic boundary 
conditions are employed, the transmission coefficients are obtained after their integration 
over the /c-points. The current / as a function of VJ, is obtained from 

/(H) = I jdE T{E, H) [f{E + eH/2) - f{E - eV,/2)], (2) 

where f{E + eVb/2) and f{E — eVb/2) are the Fermi-Dirac distribution functions of the left 
and right electrodes. The electronic temperature used is 10 K. 

III. RESULTS AND DISCUSSION 

We examine the nine interface geometries (Fig. [2]) leading to the five bonding types and 
two molecular orientations that were discussed in Sec. I. Their molecular orientations, linker 
groups, bonding types, binding sites, and distances d between the electrodes are listed in 
Table 1. 

A. Interface geometries considered 

In the Au-S bonding, four interface geometries are considered: Fig. [2t^a), (b), (e), and 
(f). In Fig. [2]^a), Geo l:Au-(SC3ll6)2-hollow, one SCsHg linker molecule is bonded to each 
Au electrode, where the S atom is adsorbed at the hollow sites of the flat Au surface. For 
a fixed distance d, after the geometry relaxation, this geometry bears the shortest distances 
between the S and Au surface atoms in the range of 2.52 to 2.82 A. If the S atom is bonded 
to the on-top site of the Au surface [Geo l:Au-(SC3H6)2-ontop, Fig. EJ^b)], the total energy 
increases by 0.094 eV compared to the hollow-site case, and the shortest distance between 
the S and Au atoms becomes 1.92 A. In Fig.El^e) and (f), two S linker atoms [Geo 2:Au-S4] 
or two SC2H4 linker molecules [Geo 2 :Au-(S 02114)4] are bonded to the hollow sites of the 
surface of each electrode. For these two geometries the shortest distances between the S and 
Au atoms are in the range of 2.52 to 2.68 A. 

Our study of the Au-C bonding was motivated by recent measurements on the con- 
ductance of a molecular junction based on benzene directly bonded to Pt electrodes.— It 
is prevalent that chemical bonding between a metallic surface and C atoms is difficult to 
achieve under ordinary conditions, especially between an Au surface and C atoms. However, 
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the direct Au-C bonding does occur in a few exceptional cases, such as isocyanide derivatives 
adsorbed on Au surfaces.—"^ In the geometry with the Au-C bonding [Geo l:Au-(C3H6)2, 
Fig. [21(c)], one CsHg hnker molecule is adsorbed at the on-top site of each electrode, and the 
shortest distance between the Au and outermost C atoms is 2.28 (2.37) A for the left (right) 
linker molecule. The study of the Au-Au bonding was inspired by an experimental effort 
to increase conductance through small molecules using Au-Au bonding.— In the geometry 
with the Au-Au bonding [Geo l:Au-(AuC3H6)2, Fig. Et^d)], one AuCaHg linker molecule is 
bonded to the hollow sites of each electrode, and the shortest distances between the Au 
surface and the Au atoms from the linker molecules are in the range of 2.73 to 2.79 A. 

The study of the geometries with physisorption and without linker molecules, was moti- 
vated by the transport measurement through the Mni2 without any linker molecules^. We 
consider such geometries in the orientation (2). In the geometry with the Au-H bonding 
[Geo 2: Au-H- no- linker. Fig. |2|^h)], no linker molecules are placed between the Mni2 and 
the electrodes. However, the short distance d in that geometry allows Au-H bonding to 
be initially formed at four different locations. We choose one of the four Au-H bonding 
distances in its initial geometry to be slightly shorter than the three distances, and then 
upon geometry relaxation, only one Au-H bonding remains effective as shown in Fig. EJ^h). 
In the case of physisorption, two geometries are considered: Figs. [2](g) and (i). In Fig. [2|^g), 
Geo 2:phys-(CH3)4, two CH3 linker molecules are attached to each side of the Mni2. Upon 
geometry relaxation (with a fixed d until the magnitude of the maximum force is less than 
0.1 eV/A), in this geometry, the shortest distance between the Au and the C atoms from 
the linker molecules is 2.41 A, which is long enough such that chemical bond is not formed 
between the linker molecules and the electrodes. In Fig. [2]^i), Geo 2:phys-no-linker, no linker 
molecules are placed between the Mni2 and the electrodes, and the Au surface atoms are 
well separated from the H atoms from the Mni2 in the range of 2.97 to 3.10 A. Thus, no 
chemical bond is formed between the electrodes and the Mni2 in this geometry, either. The 
geometry of Geo 2 :phys- no-linker was built from the optimized geometry of Geo 2:Au-H-no- 
linker with an increased distance d. Hereafter, when a particular bonding type is discussed, 
its corresponding binding site is also implicitly assumed. 
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FIG. 3: Majority-spin density of states projected onto the Mn d orbitals of the (a) initial and (b) 
optimized geometries, Geo l:Au-(SC3He)2-hollow. Refer to Fig. [2] for numbering of the Mn sites. 
The densities of states for Mn(3) and Mn(4) sites are the same as those for Mn(l) and Mn(2). 

B. Effect of geometry relaxation 



For an isolated Mni2 molecule in the gas phase, without the on-site Coulomb repulsion U, 
our DFT calculations using SIESTA show that the energy gap between the highest occupied 
molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) levels is 
0.38 eV. Since the Mni2 has a high magnetic moment of 20 fiB in the ground state, both the 
HOMO and LUMO arise from the majority-spin Mn d orbitals. This feature of the HOMO 
and LUMO does not change with inclusion of a proper value of the U term, although the 
HOMO-LUMO energy gap increases with the U term. This was shown in previous DFT+U 
calculations^!^, in contrast to the result discussed in Ref.[l6|. Such a characteristic of the 
HOMO and LUMO was also earlier revealed by all-electron DFT calculations^, and it was 
indirectly demonstrated in experiments on locally charged Mni2 molecules^*^. 

When an Mni2 is bridged between the electrodes in the orientation (1), the S4 symmetry 
of an isolated Mni2 is broken and the molecular levels in the interface show approximately 
twofold symmetry. In addition, the magnetic easy axis of the Mni2 is slightly tilted from the 
axis perpendicular [x axis. Fig. [2t^a)] to the transport direction, in order to allow the linker 
molecules to be bonded to the electrodes. Furthermore, geometry relaxation at a fixed d 
renders some noticeable changes in the molecular levels. To compare the molecular levels 
before and after the geometry relaxation, we examine the geometry Geo l:Au-(SC3H6)2- 
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FIG. 4: (Color online) (a) Majority-spin transmission T{E,0) at zero bias as a function of energy 
E relative to the Fermi level Ej for the five interface geometries with the orientation (1), and 
(b) for the five geometries with the orientation (2). The arrows indicate the transmission peaks 
associated with the HOMO and LUMO levels of the Mni2. 

hollow. For its initial geometry, the HOMO and LUMO levels of the Mn d orbitals, are 
separated by 0.3 eV, which is similar to the HOMO-LUMO energy gap of 0.38 eV for the 
isolated Mni2. With the geometry optimization, the Mn d orbital levels right below the Fermi 
level Ef approach toward Ef, but the Mn d levels right above Ef become split (Figj3]). As 
a result, the separation between the HOMO and LUMO levels reduces to 0.15 eV. The 
HOMO level is located near —0.15 eV, while the LUMO level is near Ef [Fig. [3|b)] (this is 
due to some charge transfer from the Au electrodes to the Mni2^). The same kind of shift 
is expected for the rest of the geometries considered in the orientation (1). Now when an 
Mni2 is bridged between the electrodes in the orientation (2), the S4 symmetry is by large 
preserved in the molecular levels in the interface. If the distance d is not too short, the 
geometry relaxation will shift the molecular levels for geometries in the orientation (2) as 
well. One exception is Geo 2:Au-S4, where a combination of the strong Au-S bonding with 
the short distance d prevents the Mni2 and the four S atoms from being significantly rotated 
or stretched/compressed with geometry relaxation. 

We discuss the effect of geometry relaxation on T{E,0). Similarly to the effect on the 
molecular levels, the geometry relaxation gives rise to a upward shift in the transmission 
peaks right below Ef and more spread peaks right above Ef, as shown in the first and the 
third panels of Fig. HJ^a) from the top. According to our calculations, the current is carried 
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via the LUMO level of the Mni2, and thus the broadening of the LUMO level is relevant to 
the current-voltage characteristics. Without contributions from interactions with phonons, 
in resonance, the individual transmission peak centered at Ei for given VJ, (for a single 
conductance channel) is solely described by the broadening of the corresponding molecular 
level as follows:—'^ 

where Ti{Ei,Vb) and T2{Ei,Vb) are the broadening of the ith molecular level caused by 
the left and right electrodes, respectively, and they are determined by full widths at half 
maximum of the transmission peak centered at Ei. Here T{Ei, Vb) is the average broadening 
given by (Ti{Ei, Vb) + T2{Ei, Vb))/2. (Notice that for symmetric coupling, Fi = Eq. ([3]) 
implies that T becomes unity at E = E^.) The orbital broadening Fi 2(-E'i, Vb) is related to 
the hopping integral ti^2 (between the electrodes and the Mni2) as 7itl 2P{Ei), where p{Ei) 
is the density of states at the ith level.— The values of ti,2 and p{Ei) depend on bonding 
types and binding sites. We compute the values of Fi and F2 from fitting the transmission 
peaks to Eq. ([3]). For Geo l:Au-(SC3He)2-hollow, the geometry relaxation reduces the value 
of Flumo from 0.0028 to 0.0020 eV, while it enhances the value of Fhomo from 2.6x10"^ to 
5.0x10"^ eV, as listed in Table 1. Since the value of Flumo decreases to a small degree with 
the geometry relaxation, a slightly lower current flows through the relaxed junction than via 
the initial one (the topmost left panel of Fig. [5]). Similar behavior is expected for the rest of 
the geometries. Henceforth, we discuss optimized interface geometries only, unless specified 
otherwise. 



C. Effect of molecular orientation 



For the interface geometries with the orientation (1), Fig. [2]^a), (b), (c), and (d), the 
linker molecules are closer to the Mn(5), Mn(6), Mn(9), and Mn(lO) sites than the rest of 
the eight Mn sites, as illustrated in Fig. El Among the former four Mn sites, the Mn(5) 
and Mn(9) sites predominantly contribute to the LUMO of the Mni2 [Fig. EJ^a)], as well as 
to some degree to the HOMO [Figs. (Kb) and[6](b)]. Notice that the HOMO and LUMO 
arise from the majority-spin Mn d orbitals. A much greater contribution to the HOMO 
arises from the Mn(8) and Mn(12) sites [Fig. Et^b)], but these sites are not as close to the 
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FIG. 5: (Color online) Current-voltage (I—V) characteristics for six interface geometries: majority- 
spin contribution (red box), minority-spin contribution (blue star). The topmost left panel shows 
I — V curves for both the initial (brown circle) and the optimized geometries Geo l:Au-(SC3H6)2- 
hollow (red box). Note that Geo l:Au-(SC3He)2-hollow and Geo 2:phys-(CH3)4 have different 
vertical scales from the rest of the geometries. 

linker molecules as the former four sites. Thus, for the geometries with the orientation 
(1), the broadening of the LUMO level is expected to be somewhat larger than that of the 
HOMO level, as shown in Table 1 and Figs. Et^b), HJ^a), and [71 For the geometries with 
the orientation (2), in the case of chemisorption. Fig. |2]^e), (f), and (h), the linker molecules 
are in a closer proximity to the Mn(l), Mn(2), Mn(3), Mn(4), Mn(6), Mn(8), Mn(lO), and 
Mn(12) sites, than to the other four Mn sites, as shown in Fig. [2l The Mn(6), Mn(8), 
Mn(lO), and Mn(12) sites predominantly contribute to the HOMO [Fig. IH^d)] rather than 
the LUMO [Fig. Et^c)], and so the broadening of the HOMO level is much larger than that 
of the LUMO level [Table 1 and Fig. Hl^b)]. In Fig. lU^b), the transmission peak associated 
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FIG. 6: Spatially resolved density of states (DOS) projected onto (a) the LUMO and (b) HOMO 
for Geo l:Au-(C3H6)2, and onto (c) the LUMO and (d) HOMO for Geo 2:Au-H-no-linker. The 
LUMO and HOMO originate from the majority-spin Mn d orbitals. The criterion of the isosurface 
shown is 2.0 x 10~^ electrons/a^, where is the Bohr radius. The numbering of the Mn ions is 
based on Fig. [2l The blobs far away from the core of the Mni2 (marked by L) in (a) and (b) arise 
from the linker molecules. 

with the HOMO level for Geo 2:Au-S4 (Geo 2:Au-H-no-linker), is very broad because it 
corresponds to four (two) molecular levels broadened right below Ef including the HOMO 
level. In the case of physisorption, Figs.|2]^g) and (i), the Mn(6), Mn(8), Mn(lO), and Mn(12) 
sites are much closer to the electrodes than the Mn(5), Mn(7), Mn(9), and Mn(ll) sites. 
The HOMO originates mainly from the former four Mn d orbitals, while the LUMO comes 
from the latter four Mn d orbitals. Thus, for the geometries in the orientation (2), either in 
the case of chemisorption or physisorption, the broadening of the HOMO level is expected 
to be much larger than that of the LUMO level [Table 1 and Fig. Ht^b)]. 

Now we compare the level broadening for the geometries in the orientation (1) to that 
in the orientation (2), for a given bonding type and fixed distance d. For this comparison, 
we examine the following three geometries with the Au-S bonding: one geometry in the 
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FIG. 7: Majority-spin density of states (DOS) projected onto Mn d orbitals for (a) Geo l:Au- 
(C3H6)2 and (b) Geo l:Au-(AuC3H6)2- Refer to Fig. [2] for numbering of the Mn ions. The 
densities of states for Mn(3) and Mn(4) sites are the same as those for Mn(l) and Mn(2) sites. 

orientation (1), Geo l:Au-(SC3H6)2-hollow, and two geometries in the orientation (2), Geo 
2:Au-S4 and Geo 2:Au-(SC2H4)4. The distance d in the former geometry, 25.73 A, is much 
longer than those for the latter two geometries. Thus, we extrapolate to find the broadening 
of a corresponding geometry in the orientation (2) with (i=25.73 A. The Mni2 and the linker 
molecules play the role of a potential barrier in the electron transport, and so the broadening 
of the relevant molecular levels decays exponentially as a function of d. Thus, we apply 
rLUMO=Ci exp(— C2(i), where Ci and C2 are positive constants, to the two geometries in 
the orientation (2). Extracting the values of Ci and C2, we find that the extrapolated value 
of 

Tlumo for rf=25.73 A is 0.0015 eV, which is less than the actual value of Flumo for 
the first geometry in the orientation (1), 0.0020 eV. When the same logic is applied to the 
HOMO level, the extrapolated value of Fhomo for such a geometry becomes 0.0069 eV,— 
which is much greater than the actual value of Fhomo for the first geometry, 0.00050 eV. 
Consequently, for a given bonding type and fixed distance d, the broadening of the LUMO 
(HOMO) level for geometries with the orientation (1) is more (less) pronounced than that 
with the orientation (2). Combining this result with the finding that the LUMO level 
broadens more than the HOMO level for the geometries with the orientation (1) and the 
opposite holds for those with the orientation (2), we reach the following conclusion. Although 
more statistics is desirable and other bonding types can be tested, our calculations indicate 
that as long as the current flows through the LUMO level, (for a given d and bonding type) 
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FIG. 8: (Color online) Spin-polarized DOS projected onto the Mn d orbitals for (a) Geo l:Au- 
(SC3H6)2-hollow and (b) Geo l:Au-(C3He)2: majority-spin (black) and minority-spin (red). Refer 
to Fig. [2] for numbering of the Mn ions. The DOS projected onto the Mn(l) and Mn(2) sites are 
similar to those onto the Mn(3) and Mn(4) sites. 

geometries with the orientation (1) provide somewhat a higher current than those with the 
orientation (2). However, if one can arrange the transport to occur through the HOMO 
level (by application of gate voltage to the Mni2), geometries in the orientation (2) will give 
a much higher current than those in the orientation (1). One caveat is that for a short 
distance d the size of the Mni2 prohibits interface geometries with the orientation (1) from 
being formed. 

For all of the geometries considered, the effect of spin filtering remains robust. The spin- 
filtering effect occurs in the transport through the Mni2 because its majority-spin HOMO 
and LUMO levels are well separated from the minority-spin HOMO and LUMO levels and 
because the minority- spin HOMO level is located sufficiently below the majority-spin HOMO 
level (Fig. 2 in Ref.[15|, Figs. El El and [10]) . This feature of the HOMO and LUMO does 
not change with interface geometry. However, the degree of the spin filtering depends on 
molecular orientations, bonding types, and the distance d. We first compare the spin-filtering 
effect for the geometries in the orientation (1) to that in the orientation (2). Figure [TT] shows 
the spin-polarized T{E,0) for six geometries [Fig. Et^a), (c), (d), (e), (h), and (i)]. In the 
last column of Table 1, the ratio of the majority-spin to the minority-spin transmission 
coefficients at the energy level associated with the LUMO level, {T^ — T^/T^, is provided. 
For the geometries in the orientation (1), the ratio varies from 6.11x10^ to 5.26x10^, while 
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FIG. 9: (Color online) Spin-polarized DOS projected onto the Mn d orbitals for (a) Geo 2:Au-S'4 
and (b) Geo 2:Au-H-no-linker: majority-spin (black) and minority-spin (red). Refer to Fig. [2] for 
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FIG. 10: (Golor online) Spin-polarized DOS projected onto the Mn d orbitals for Geo 2:Au- 
(302114)4: majority-spin (black) and minority-spin (red). Refer to Fig. [2]for numbering of the Mn 
ions. The DOS projected onto the Mn(l) and Mn(2) sites are similar to those onto the Mn(3) and 
Mn(4) sites. 

for those in the orientation (2), the ratio is in the range of 7.39x10^ to 2.86x10^. A larger 
spin-filtering effect is obtained for the geometries in the orientation (1) than for those in the 
orientation (2). This is partly because the linker molecules in the orientation (2) are in a 
closer proximity to the Mn(l), Mn(2), Mn(3), and Mn(4) sites than those in the orientation 



19 



(1). These four Mn sites mainly contribute to the minority-spin density of states as shown 
in Figs. E] and [91 Due to this proximity of the four Mn sites, the minority-spin DOS for 
the geometries in the orientation (2) is more spread than that for the geometries in the 
orientaion (1). An additional reason is as follows. As discussed earlier and shown in Figs. El 
[HI anddni the LUMO arises from the majority-spin Mn d orbitals. At the energy level where 
the ratio is obtained, the majority-spin transmission is of the order of unity, independent of 
the distance because the transport is carried by the LUMO level. However, at that energy 
level, the transport of minority-spin electrons is in the tunneling regime, since there are no 
corresponding molecular levels of the minority-spin. Thus, the minority-spin transmission 
decreases exponentially with d. Since the geometries in the orientation (1) have much 
longer distances d than those in the orientation (2), the coefficients of the minority-spin 
transmission for the former geometries are much smaller than those for the latter geometries. 
This renders the greatly enhanced ratio for the geometries in the orientation (1). Now let us 
briefly discuss the influence of bonding type on the spin-filtering effect. For the orientation 
(1) the minority-spin DOS for the Au-C bonding is slightly broader than that for the Au-S 
bonding [compare Fig. [Hl^a) to (b)]. Thus, the ratio for the Au-C bonding is lower than that 
for the Au-S bonding (Table 1) for a given distance d. Similarly, for the orientation (2), 
the minority-spin DOS for the Au-S bonding is more delocalized than that for the Au-H 
bonding [compare Fig. [9t^a) to (b)]. So the ratio for the Au-S bonding is lower than that for 
the Au-H bonding for a given distance d. For a specific bonding type, the ratio increases 
with increasing the distance d because of the reason explained above. For example, for the 
Au-S bonding, the minority-spin PDOS for Geo 2:Au-(SC2H4)4 is much more localized than 
that for Geo 2:Au-S4 [compare Fig. [T0]to|9]^a)]. Thus, the ratio for the former is two orders 
of magnitude greater than that for the latter. Additionally, we note from Figs. HI and [TT] that 
this ratio, {T^ — T^)/T^, remains in the range of 100 to 10^ even at the energy level where 
the transport of both the majority spin and the minority spin is in the tunneling regime. 
The spin-filtering effect may be achieved for other SMMs with high magnetic moments and 
stable ground-state spin multiplets. 
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FIG. 11: (Color online) Spin-polarized transmission at zero bias with a logarithmic scale in the 
vertical axis. Majority-spin contribution: red, minority-spin contribution: blue. The vertical scales 
of the figures on the left differ from those on the right. 

D. Effect of bonding type and linker group 



We investigate the effect of binding site on transport properties for a given bonding type, 
for example, the Au-S bonding in the orientation (1), Geo l:Au-(SC3H6)2- We consider 
hollow and on-top binding sites of the S atoms to the Au surface in that geometry, with- 
out geometry relaxation. The broadening of the LUMO level somewhat decreases from the 
hollow-site to the on-top geometries, while the broadening of the HOMO level slightly in- 
creases (Table 1). This difference in the broadening is much smaller than that induced by 
different bonding types. Despite this small difference, the overall features of T{E, 0) for the 
two geometries are fairly similar to each other, as shown in Fig. Hl^a). Thus, the I-V curve 
for the on-top site geometry is expected to be similar to that for the hollow-site geometry, 
Geo l:Au-(SC3H6)2-hollow, shown in Fig. |5l 

We compare the transport properties of the geometries with three bonding types in the 
orientation (1), where the distance d is fairly similar to one another: Au-C bonding [Geo 
l:Au-(C3H6)2], Au-S bonding [Geo l:Au-(SC3H6)2-hollow], and Au-Au bonding [Geo l:Au- 
(AuC3H6)2]- The broadening of the LUMO level for the geometry with the Au-C bonding 
is one order of magnitude larger than that with the Au-S bonding, and it is twice as large 
as that with the Au-Au bonding, as shown in Figs. HJ^a) and [71 Accordingly, the current 
for the Au-C bonding at 50 mV is one order of magnitude higher than that for the Au-S 
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bonding, and is twice as high as that for the Au-Au bonding (Fig. |5]). Among all of the 
interface geometries considered, for a fixed distance ci, the geometry with the Au-C bonding 
reveals the largest broadening of the LUMO level, leading to the largest current (Fig. 

To understand the nature of the bonding between the linker molecules and the electrodes, 
we compute the density of states projected onto the p orbitals of the S or C atoms of the 
linker molecules closest to the electrodes, and onto the s orbitals of the three Au surface 
atoms closest or bonded to those S or C atoms. The DOS projected onto the Au s orbitals 
change with interface geometries due to the interactions between the linker molecules and 
the electrodes. For Geo l:Au-(C3H6)2, only the Pz orbitals of the C atoms of the linker group 
greatly overlap with the Au s orbitals of the electrodes near Ef [the bottom four panels of 
Fig. [T^ a)]. The density of the LUMO coincides with the peaks in the region bound by the 
dashed lines in Fig. [T2]^a) . Strong overlap (in terms of the peak height and the width of the 
region) among the Au s orbitals, the C p orbitals, and the LUMO results in a high current 
through the Mni2. For Geo l:Au-(SC3H6)2-hollow, the Px, Py, and Pz orbitals of the S atoms 
have much weaker overlap with the Au s orbitals in the region where the LUMO appears 
(confined by the dashed lines) than the Au-C bonding case [the top four panels of Fig. [T2l^a)] , 
which gives rise to a reduced transmission probability. For Geo 2:Au-S4, only the p^ and Py 
orbitals of the S atoms bear some overlap with the Au s orbitals near Ef (in the region bound 
by the dashed lines) [the top four panels of Fig. [TW b)]. This overlap is weaker than that for 
the Au-C bonding but stronger than that for Geo l:Au-(SC3H6)2-hollow or for Geo 2:Au- 
(SC2H4)4 [the bottom four panels of Fig. [T2l(b)]. This tendency agrees well with the trend in 
the value of Flumo and the current. Additionally, our calculated spatially resolved density 
of states of the LUMO [Fig. EJ^a)] uncovers that the geometry with the Au-C bonding gives 
rise to larger density in the linker molecules than any other geometries considered. Thus, the 
transmission probability is highest for the geometry with the Au-C bonding. The geometries 
with the Au-S bonding do not have as high transmission probabilities as that for the Au- 
C bonding or the Au-Au bonding. Consequently, the Au-C bonding provides the highest 
current and the Au-S bonding gives rise to a lower current than the Au-C bonding or the 
Au-Au bonding. 

We discuss the transport properties of the geometries with physisorption and without 
linker molecules. One of the geometries without linker molecules, Geo 2 :Au-H- no-linker, has 
the shortest distance d among the geometries considered. The broadening of the LUMO level 
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FIG. 12: Majority -spin DOS projected onto the s orbitals of tlie three Au surface atoms and px, Py, 
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and Geo l:Au-(C3He)2, and (b) Geo 2:Au-S4 and Geo 2:Au-(SC2H4)4. In the region confined by 
the dashed lines, the density peaks overlap with the LUMO. The vertical scale in this figure differs 
from those in Figs. [3] and [71 

for that geometry is 54% of that with the shortest Au-S bonding (Geo 2:Au-S4), and it is 21% 
of that with the Au-C bonding, as shown in Table 1. The current for the geometry with the 
Au-H bonding at 50 mV is 44% of that for Geo 2:Au-S4, and it is 36% of that with the Au-C 
bonding (Fig. [5]). In the case of physisorption, with the fixed distance the broadening of 
the LUMO level for Geo 2:phys-no-linker is the same as that for Geo 2:phys-(CH3)4. This 
implies that the alkane chains do not play a role in the broadening of the LUMO level. For 
the given distance ci, the broadening of the LUMO level for these two physisorbed geometries 
takes up only 15% of that for the Au-S bonding (Geo 2:Au-S4), as listed in Table 1. Yet, the 
current for Geo 2:phys-(CH3)4 at 50 mV is two orders of magnitude lower than that for Geo 
2:Au-S4 (Fig. [5]). Thus, in this case, the ratio of the values of Flumo for the two geometries 
does not have the same order as the ratio of the corresponding currents. The discrepancy 
between the two ratios arises from the difference in the response of the junctions or in the 
changes of the transmission with bias voltage. Notice that Fig. H] shows the transmission 
spectra for zero bias. For chemisorbed junctions the molecular orbitals strongly overlap with 
the orbitals of the electrodes. Thus, when bias voltage is applied, this overlap is enhanced 
especially in the regions between adjacent transmission peaks. However, for physisorbed 
junctions, the transmission peaks remain sharp even with bias voltage because of no overlap 
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of the molecular orbitals with the orbitals of the electrodes. A combination of this effect 
with the ratio of the values of Flumo explains the ratio of the current through Geo 2:phys- 
(CH3)4 to that through Geo 2:Au-S4. Our overall comparison of the transport behavior 
among different bonding types summarizes that the current through the Mni2 decreases in 
the following order, for a given distance d: Au-C, Au-Au, Au-S, Au-H, and physisorption. 



IV. CONCLUSION 



We simulated a single-molecule junction based on the SMM Mni2 connected to Au(lll) 
electrodes via five different bonding mechanisms and two molecular orientations, and cal- 
culated their transport properties, using the non-equilibrium Green's function method and 
the spin-polarized DFT. Geometry relaxation of the junction renders a small decrease in 
the electric current at low bias voltage. If the electron transport is carried by the LUMO 
(HOMO) level of the Mni2, geometries in the orientation (1) provide a higher (lower) current 
that those in the orientation (2). This is due to the structure of the Mni2. Our two-terminal 
transport calculations reveal that the LUMO level is relevant to the transport through the 
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n all of the interface geometries considered, the spin-filtering effect discussed in 
Ref. jl4| remains robust, and it may occur to some other SMMs where the ground-state 
spin is large and the ground-state spin multiplet is reasonably well separated in energy from 
the low-lying excited spin multiplets. Some experimental studies^i^i^>ii show that upon de- 
position of Mni2 molecules on an Au surface, the valence states of all of the Mn ions are 
preserved, while some other experimental works^i^ii^ reveal that such a deposition induces 
changes in the valence states of some of the Mn ions. In cases where such changes do not 
occur-i^ii^iii, the spin-filtering effect is expected. We also find that the current through the 
Mni2 depends on a bonding type and that it decreases in the following order for a fixed 
d: Au-C, Au-Au, Au-S, Au-H, and physisorption. This is because the overlap among the 
LUMO, the S or C p orbitals from the linker molecules, and the s orbitals of the Au surface 
atoms, decreases in that order. The degeneracy in the magnetic levels of the Mni2 for a 
given spin S is lifted by spin-orbit coupling. Then inelastic transport through low-energy 
spin excitations must be included for quantitative comparison to experimental data. How- 
ever, the general trend in the transport as a function of bonding type and interface geometry 
found in this study may still hold. 
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TABLE I: Average broadening of the HOMO and LUMO levels (rnoMO, Tlumo) of the Mni2 
for the nine interface geometries. All geometries are optimized, unless stated otherwise. The 
broadening is calculated from zero-bias transmission spectra. The value of Fhomo for Geo 2:Au-S4 
(Geo 2:Au-H-no-linker) is marked by f because it represents the full width at half maximum of 
the single broad transmission peak associated with the four (two) molecular levels right below Ef 
including the HOMO level, as shown in Fig. \Mi>)- The ratio of the majority-spin to the minority- 
spin transmission, (T-[- — T^j/T^, is computed at the energy level corresponding to the LUMO 
level. 



interface geometry 
(orientation, linker group) 


bonding type 
binding site 


d{k) 


Fig. 1 


Tlumo 
(eV) 


Thomo 
(eV) 


(Tt - 


Geo l:Au-(SC3H6)2-hollow 
initial geometry 


Au-S, hollow 


25.73 


N/A 


0.0028 


2.6x10"^ 


2.57x10^ 


Geo l:Au-(SC3He)2-ontop 
initial geometry 


Au-S, on-top 


25.73 


(b) 


0.0020 


3.4x10-5 


5.26x10^ 


Geo l:Au-(SC3H6)2-liollow 


Au-S, hollow 


25.73 


(a) 


0.0020 


0.00050 


4.37x10^ 


Geo l:Au-(C3H6)2 


Au-C, on-top 


23.33 


(c) 


0.020 


0.020 


6.11x10^ 


Geo l:Au-(AuC3H6)2 


Au-Au, hollow 


25.73 


(d) 


0.0090 


0.0028 


1.40x10^ 


Geo 2:Au-S4 


Au-S, hollow 


14.48 


(e) 


0.0078 


o.iit 


7.39x10^ 


Geo 2:Au-(SC2H4)4 


Au-S, hollow 


19.01 


(f) 


0.0040 


0.016 


2.86x10^ 


Geo 2:phys-(CH3)4 


physisorption 


14.56 


(g) 


0.0010 


0.0032 


2.32x10^ 


Geo 2:Au-H-no-linker 


Au-H, on-top 


12.69 


(h) 


0.0042 


0.054t 


4.63x10^ 


Geo 2:phys-no-linker 


physisorption 


14.48 


(i) 


0.0011 


0.0092 


2.74x10^ 
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